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It is shown that different couples of stress-energy and spin tensors of quantum relativistic fields, 
which would be otherwise equivalent, are in fact inequivalent if the second law of thermodynamics 
is taken into account. The proof of the inequivalence is based on the analysis of a macroscopic 
system at full thermodynamical equilibrium with a macroscopic total angular momentum and a 
specific instance is given for the free Dirac field, for which we show that the canonical and Belin- 
fante stress-energy tensors are not equivalent. For this particular case, we show that the difference 
between the predicted angular momentum densities for a rotating system at full thermodynamical 
equilibrium is a quantum effect, persisting in the non-relativistic limit, corresponding to a polar- 
ization of particles of the order of huj/KT {uj being the angular velocity) and could in principle 
be measured experimentally. This result implies that specific stress-energy and spin tensors are 
physically meaningful even in the absence of gravitational coupling and raises the issue of finding 
■ the thermodynamically right (or the right class of) tensors. We argue that the maximization of the 

thermodynamic potential theoretically allows to discriminate between two different couples, yet for 
(~| , the present we are unable to provide a theoretical method to single out the "best" couple of tensors 

r-{ . in a given quantum field theory. The existence of a non- vanishing spin tensor would have major 

consequences in hydrodynamics, gravity and cosmology. 

^ ■ 

I. INTRODUCTION 

It is commonly known that stress-energy and spin tensors are not uniquely defined in field theory as long as 
^ • gravity is disregarded. In quantum field theory, distinct stress-energy tensors differing by the divergence of a rank 3 
1^ ' tensor provide, once integrated in three-dimensional space, the same generators of space-time translations provided 
. that the flux of the additional rank 3 tensor field (hereafter referred to as superpotential) vanishes at the boundary. 
Correspondingly, in classical field theory, the spatial integrals of T^^ yield the same values of total energy and 
^■f-^ momentum. Within field theory on a flat spacetime, it is then possible to generate apparently equivalent stress- 
^ ■ energy tensors which are e.g. symmetric or non-symmetric. Indeed, gravitational coupling provides an unambiguous 
^-H ] way of defining the stress-energy tensor; in General Relativity, it is symmetric by construction and the spin tensor 
'/"J vanishes. However, in a likely extension known as Einstein-Cartan theory (not excluded by present observations) the 
\ spin tensor is non-vanishing and the stress-energy tensor is non-symmetric. 

Can we say something more? In classical physics we have a stronger requirement with respect to a quantum theory: 
we would like the energy, momentum and angular momentum content of any arbitrary macroscopic spatial region 
to be well defined concepts; otherwise stated, we would like to have objective values for the energy, momentum and 
angular momentum densities. If these quantities are to be the components of the stress-energy and spin tensors, such 
a requirement strongly limits the freedom to change these tensors. It is crucial to emphasize, from the very beginning, 
the difference between quantum and classical tensors. The quantum stress-energy and spin tensors, henceforth denoted 
with a hat ^, are operatorial expressions depending on the microscopic quantum field operators ^I', whereas the classical 
ones are c- numbers. The relation between them is 



X 



T^"'{x)=tT[p ■.f^'^x)-] S^'^'^ix) =It[p ■.S^'^"'{x)■] (1) 

where p is the density operator describing the (mixed or pure) quantum state and : denotes normal ordering; the 
latter is usually introduced in the mean value definition in order to subtract the zero-point infinities ^. According 
to ([1} , a change of quantum stress-energy and spin tensors could induce a change of the corresponding classical ones 
in an undesirable fashion, meaning that energy or momentum or angular momentum density get changed. However, 
the change that classical mean values undergo as a reflection of a variation of quantum tensors crucially depends on 
the physical state p. Particularly, we will see that the freedom of varying the stress-energy and spin tensors at a 
quantum level depends on the symmetry features of the physical state: a highly symmetric state allows more changes 
of quantum tensors than a state with little symmetry does. 



^ For a discussion of the meaning of normal ordering for interacting fields see e.g. ref. Q|. We stress that the results obtained in this work, 
particularly in Sect. IVII are anyhow independent of the use of normal ordering in eq. JTJ. 
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In this paper, we prove that a system at full thermodynamical equilibrium with a macroscopic value of angular 
momentum, thence rigidly rotating allows to discriminate between different quantum spin tensors, and, conse- 
quently, between different quantum stress-energy tensors. This kind of inequivalence shows up only for a rotating 
system whereas all quantum tensors are equivalent for a system at the more familiar thermodynamical equilibrium 
with vanishing macroscopic angular momentum. The paper is organized as follows: in Sect. we will discuss the 
general class of stress-energy tensor transformations ensuring the invariance of conservation equations; in Sect. IIIII we 
will discuss the usual thermodynamical equilibrium distribution and its symmetries and in Sect. IIVI we will do the 
same for a system at full thermodynamical equilibrium with angular momentum; in Sect. |V]we will obtain the most 
general form of mean stress-energy and spin tensor for a system at full thermodynamical equilibrium with angular 
momentum and show that the equivalence between different quantum tensors no longer applies unless peculiar con- 
ditions are met; in Sect. IVII we will present and prove a concrete instance of inequivalence for the free Dirac field; 
finally, in Sect. IVlIl we will summarize and further illustrate the obtained result and discuss the possible consequences 
thereof. 



Notation 

In this paper we adopt the natural units, with h = c = K = 1. 
The Minkowskian metric tensor is diag(l, —1, —1, —1); for the Levi-Civita symbol we use the convention e"^^'^ = 1. 
We will use the relativistic notation with repeated greek indices assumed to be saturated. Operators in Hilbert space 
will be denoted by an upper hat, e.g. R, with the exception of the Dirac field operator which is denoted with a capital 
^. 

II. TRANSFORMATIONS OF STRESS-ENERGY AND SPIN TENSORS 

The conservation equations ensuing from the relativistic translational and Lorentz invariance are the well known 
continuity equations of energy-momentum and total angular momentum: 

df.f'"' = 

dxJ^^'"' = dx (S^^^'" + x'^f^" - x'^f = dxS^^'"' + f'"' - f '•'^ = (2) 

However, stress-energy T and spin tensor S are not uniquely defined in quantum field theory; once a particular couple 
(T, S) of these tensors is found, e.g. applying Noether's theorem to some Lagrangian density (the so-called canonical 
tensors), it is possible to generate new couples (T',5') through the following pseudo-gauge transformation [3|: 

f'fJ-f — J'/Ji' —da (^"•'^'^ — ^f^^°'^ — ^'^•"f^^ 

where $ is an arbitrary tensor of rank three antisymmetric in the last two indices depending on the fields ip. It is 
easy to check that the new couple fulfills the same continuity equations ^ as the original one and that the new total 
angular momentum tensor, like the stress-energy tensor, differs from the original one by a divergence: 



The spatial integrals over the domain D,: 



P" = / d^x T"'' 
Jn 



= / d^x J°'^'•' = / d^x S^"-' + x^f^" - x^f"" (5) 



are conserved (and are generators of translations and Lorentz transformations if the domain is the whole space) 
provided that the fluxes at the boundary vanish: 

dS* f'^'m = 

an 

ds + x^r" - x^r")/!, = (6) 

an 
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rii being the surface normal versor. The eqs. ([6]) are usuahy ensured by enforcing special boundary conditions for 
the fields; for instance, the familiar periodic boundary conditions for a box, or requiring some expression involving 
the field or its normal derivatives to vanish at the boundary (this will be the case in this paper, see Sect. IVip . It 
is very important to stress here that one can obtain conserved stress-energy and angular momentum tensors and a 
corresponding set of conserved generators even within a finite or bounded region, thus breaking global translational 
and rotational symmetry, provided that suitable boundary conditions are imposed on the fields; this point will be 
further discussed in Sect. I VI I 

The spatial integrals ([S]) of the new tensors are invariant, thus yielding the same generators, if the tensor $ is such 
that the following boundary integrals vanish, according to eq. ([3]) and ([4]): 

/ dS l^'^^" - SO'^" - S'^'*") rii = 

JdQ ^ ' 

J dS (d''-"" - - ^" '"^ - x" - - n, = (7) 

In fact, if the above conditions are met, the fiux integrals ^ of the new primed tensors vanish because the new 
tensors also fulfill the continuity equations. In conclusion, a pseudo-gauge transformation like ([3]) is always possible 
provided that the boundary conditions ([7]) are ensured; in this case the couple {T,S) and {T',S') are regarded as 
equivalent in quantum field theory because they give the same total energy, momentum and angular momentum, in 
the operatorial sense. 

The classical counterpart of transformation ([3]) can be calculated by applying eq. ([1]) to both sides and this obviously 
leads to: 

If the system is macroscopic, we would like the mean values of those tensors to be invariant under a transformation 
([8]), and not just their integrals. This is because energy, momentum and total angular momentum densities classically 
must take on objective values, independent of the particular quantum tensors. A minimal requirement would be the 
invariance of the aforementioned densities, that is: 

However, this is a frame-dependent requirement; a Lorenzt-boosted frame would measure a different energy-momentum 
density if only the first row of the stress-energy tensor was invariant under transformation ([8]) in one particular frame. 
We are thus to enforce a stricter requirement, namely: 

whereas, for the rank 3 angular momentum tensor, we can make a looser request: 

jfX,t.^ ^ j\.p.v _^ ^Xt^j^u _ ^A^^M (10) 

where if is a vector field. Indeed, if we limit ourselves to spatial indices /i, = 1, 2, 3, the above equation is enough 
to ensure that the angular momentum densities, with A = 0, are the same in any inertial frame. Comparing eq. ([9]) 
with eq. ([5]), we get: 

9„ ($"^^'' - ~ $'^^"'^) = (11) 

while comparing eq. ()10|) with the mean of eq. Q and taking ()lip into account, we obtain a simple condition for the 
superpotential to meet: 

(<i)"'^'' - <I>^^"'' - ^i'."-^) - x"" ($"'^^ - (J)-^'"^ - $*'^"^)] = q^^'K" - g^^Kf^ 
^ 1 (^^'^^'^ _ ^^.M"^ _ ^''.A'A _ ^i^Ap. ^ ^A,^^ ^ ^A'^^A) ^ = g^'^K" - g^^'K^ (12) 

Plugging this last result back into eq. ((TT]) one obtains: 

2da{K''g'"' - K-^g^"") = 28" Rt" - 2g'"'d ■K = 
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Contracting the indices /i and v we obtain at once that the divergence of the vector field K vanishes and so, because 
of the above equation: 

d^'K^ = (13) 

Therefore the eqs. (fTTj) and ([12]) imply that the vector field is a constant field. The possible directions of this field 
will be dictated by the symmetry properties of the system under consideration, as we will see in the next two sections. 

It should be emphasized that the conditions PT|) and do not need to apply to the quantum tensor <I>, which 
only has to meet the boundary conditions (O , as has been seen. On the other hand, if we take the mean values of ([7]) 
applying tr(p ) on both sides, the ensuing equation is a trivial consequence of the eq. (|TT|) . In fact, it may happen 
that the mean value of the superpotential $ fulfills eqs. (fTTj) and (fT2|) even though its quantum correspondent $ 
does not, because of specific features of the density operator p. In this case, the couples (T,5) and {T', S') are to be 
considered equivalent only with regard to a particular density operator, that is for a specific quantum state. 

We will see in the next two sections that the equivalence between couples of tensors crucially depends on the 
symmetry properties of the physical state p (either mixed or pure) . Particularly, we shall see that if p is the usual 
thermodynamical equilibrium operator, proportional to exp[— i7/T + iiQ/T], any quantum tensor $ will result in a 
mean value fulfilling eq. (jlip and (jl2p . This means that all possible quantum microscopic stress-energy and spin 
tensors will yield the same physics in terms of macroscopically observable quantities. 



III. THERMODYNAMICAL EQUILIBRIUM 



The familiar thermodynamical equilibrium distribution (in the thermodynamical limit V — >■ oo): 

p^^eM-H/T + ^iQ/T) (14) 

where Q denotes a conserved charge, and Z is the grand-canonical partition function 

Z = ti[exp{-H/T + fiQ/T)] 

is remarkably symmetric. It is space-time translationally invariant, since both Q and H commute with translation 
operators T(a) = exp[ia • P]. This entails that the mean value of any space-time dependent operator A{x)^ including 
stress-energy and spin tensor, are independent of the space-time position: 

tiip:A{x + a):] = tr[p : f (a)l(a;)f (a)-i :] = tr[pf (a) :l(a;) : f (a)-i] = tr[f (a)-ipf (a) :l(a;):] = tr[p : A{x) :] (15) 

where the ciclicity of the trace and the transparency of the normal ordering with respect to translations have been 
used ^. As a consequence, the mean value of any space-time derivative vanishes, and so will do the divergences on 
the right hand side of eq. Therefore, the mean stress-energy tensor will be the same regardless of the particular 
microscopic quantum tensor used. For instance, for the Dirac field, the three tensors: 

i^^^d"^ |'*7''a''* ^ ¥7^9"^' + {p^v) (16) 

will result in the same mean stress-energy tensor. 

Also, the density operator (|14p manifestly enjoys rotational symmetry, for H and Q commute with rotation operators 
R. This implies that most components of tensors vanish. To show that, it is sufficient to choose suitable rotational 
operators and repeat the same reasoning as in eq. (fT5|) . For instance, choosing the R2(7r) operator, i.e. the rotation 
of 180 degrees around the 2 (or y) axis, changing the sign of 1 (or x) and 3 (or z) components and leaving 2 and 
unchanged, in the x = (t, 0) one has: 

ri2(^) ^tr[p:fi2(x):] =tr[R2(7r)pR2(7r)-i : f ^^(a;) :] = tr[p R2(7r)-i ■.f^\x):R2{n)] 
= tr[pR2(^)iR2(^)^ :f^^(R2(7r)-i(a;)):] = -tT[p -.f^HR^inr^x)) :] ^ -iT[p -.f^^x) :] = -T^\x) (17) 



^ Here a comment is in order. The transparency of the normal ordering with respect to a conjugation transformation, that is : AF(^)A~^ := 
A :_F('J') : A~^ where A is a translation or a Lorentz transformation and F a function of the fields and its derivatives, is guaranteed for 
free fields provided that the vacuum |0) is an eigenstate of the same transformation, which is always the case. For interacting fields, we 
will assume that the definition of normal ordering (for this problem, see e.g. ref. 0) is such that transparency for conjugation holds; 
anyhow, for the examined case in Sect. IVll we will just need transparency for a free field. 
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where, in the last equahty, we have taken advantage of the homogeneity of ah mean values shown in eq. (|T5|) : thus, 
r^^((t,0)) — and, in view of the translational invariance T^'^{x) — Vx. Similarly, by choosing other rotation 
operators, it can be shown that all off-diagonal elements of a tensor vanish. The only non-vanishing components are 
the diagonal ones, which, again owing to the rotational symmetry (choose Ri(7r/2) and repeat the above reasoning), 
are equal: 

The component r''°(x) can also be non- vanishing and its value is unrelated to the other diagonal ones. Altogether, 
the mean stress-energy tensor can only have the diagonal (symmetric) form: 

/ p \ 

^"'= p =iP + P)i''i''-P9''- 

V pj 

where t is the unit time vector with components (1, 0) and p and p have the physical meaning of proper energy density 
and pressure. It should be stressed that, for a system at full thermodynamical equilibrium described by p in eq. (jl4p 
they would be the same regardless of the particular form of the quantum stress-energy tensor, e.g. those in eq. 
for the free Dirac field. 

As far as the superpotential is concerned, it is easy to convince oneself, by using suitable rotations, that the only 
non- vanishing components are: 

Hence, one scalar function B, independent of x, is sufficient to determine the spin tensor for a system at full thermo- 
dynamical equilibrium: 

This tensor has exactly the form for a "good" superpotential derived in eq. (IT^ fulfilling condition (|13l) . In conclusion, 
any transformation of the kind ([3]) will yield the same energy, momentum and angular momentum density for all 
inertial frames and so, all quantum stress-energy and spin tensors are equivalent as far as the density operator (jl4p 
is concerned. 



IV. THERMODYNAMICAL EQUILIBRIUM WITH ANGULAR MOMENTUM 

The situation is remarkably different for a thermodynamical system having a macroscopic non-vanishing total 
angular momentum. In this case, in its rest frame (defined as the one where the total momentum vanishes) the 
density operator reads [l, : 

p= ^exp{-H/T + uj-J/T + fiQ/T) (18) 

where a; has the physical meaning of a constant, fixed angular velocity around which the system rigidly rotates. The 
factor is the rotational grand-canonical partition function: 

Z„ = tr[exp(-i//r + a; • J/T + pQ/T)] (19) 

The density operator (I18p is much less symmetric than that in (I14[) and this has remarkable and interesting con- 
sequences on the allowed transformations of stress-energy and spin tensor. The surviving symmetries in (jlSI) are 
time-translations T(i) and translations along the u: axis T{z), rotations around the uj axis Rcj{(p) and reflection 11;;^ 
with respect to planes orthogonal to a; (if H is parity-invariant). 

The density operator (jlSp can be obtained in several fashions: by maximizing the entropy with the constraint of 
fixed mean value of angular momentum [S] , generalizing to the quantum-relativistic case an argument used by Landau 
for classical systems (3| or as the limiting macroscopic case of a quantum statistical system with finite volume and 
fixed angular momentum in its rest frame in an exact quantum sense, i.e. belonging to a specific representation of the 
rotation group It should be pointed out that H,Q and the angular momentum operator along the u) direction 
commute with each other, so that the exponential in (1181) also factorizes. 
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The density operator ([T8)) implies that, in its rest frame, the system is rigidly rotating with a velocity field v = u) x x. 
The classical, non-relativistic derivation by Landau |3j] shows this in a very simple fashion by assuming that the system 
is made of macroscopic cells. To show the same thing within a quantum formalism, implies a little more effort, which 
is nevertheless quite enlightening. Consider a vector field V{x) and calculate its mean value at a point .t + a by using 
space-time translation operators. Evidently: 

tT[p:V''{x + a)-] =tr[pf(a) : ^''(x) : f (0)-^] = tr[f (a)" VT(a) -.V'^ix):] 



-tr[f(a)-ie-^/^+'^-J/^+^«/^f(a) :V''{x): 



.^^[g-i?/T+a;.f(a)-ijf(a)/T+MQ/T .^-(3.).] (20) 



where known commutation properties [Q, P^] = and [H, P^] = have been used. Now, from the theory of Poincare 
group is known that: 



T(a)-ijT(a) = J + ax P 



(21) 



whence, from ([20 



tilp-.V^ix + a):] = ±tr[e-^/^+'^-(5+'^xP)/^+^Q/^ ■.V''{x):] = ±tr[e-^/^+('*'><'^)-P/^+'^-J/^^+^Q/^ -.V^ix):] (22) 
Zoj Z^ 

We can now define the temperature four-vector: 



which can also be rewritten as: 



/3 = -(l,u;xa) 



/3 = 7^" = 7^(7, 7v) 
-to -to 



(23) 



where v = a; x a, 7 = — and To = 7T. The vector v is manifestly a rigid velocity field, while Tq is the inverse 

modulus of /3, i.e. the comoving temperature which differs by the constant uniform T by a 7 factor @, 0. The mean 
value of V^(x -t- a) in eq. ([2^ becomes: 



tr[p ■V''{x + a) :] = ^tr[exp(-^(a) • P + a; • J/T + /iQ/T) : V"" {x) :] 

Zuj 

Since /3 is timelike (provided that w < 1), it is possible to find a Lorentz transformation A such that: 



(24) 



(25) 



A convenient choice is the pure Lorentz boost along the v = u) x a direction, which, being ortogonal to a;, leaves the 
operator J • a; invariant: 

A — exp[— i arccosh(7)v • K] 

where (i = 1,2,3) are the generators of pure Lorentz boosts. Thereby, the trace on the right hand side of the 
eq. (|24l) can be written: 



tr 
tr 



tr 



_X-i(pO/To+7ix.>-J/To-l-7MQ/T'o)A . [x) : 



tr 



'^-l^-P°/To+7^-J/To+7f^Q/To^ :i>''(a;) 

= (A-i)-tr[c-^°/^"+^'^-^/^«+^^'5/^° -.V^^iAix)) 
Finally, from ([24]) and ([Ml) we get: 



g-P°/To+7'^-J/To+7MQ/Toy^ .^(2;): A"^ 



1 



tr[p:V''ix + a):]^—iA-Xtr 



■.v^[h{x)Y 



(26) 



(27) 
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which tells us how to calculate the mean value of a vector field at any space-time point given its value in some other 
specific point. The most interesting feature of eq. (P?]) is that the density operator on the right hand side the same 
as p on the left hand side with the replacement: 



T — > To = 7(a)r ijj — > 7(0)0; /i — > ^{a)ii 
If we choose x — (0, 0), i.e. the origin of Minkowski coordinates, and a — (0, a) the eq. ([27]) implies: 
tr[p :i/^(0,a):] = —{\-^Yir L-^5°/t„+7-J/To+7mQ/To :i/A'(o,0):" 



(28) 



(29) 



that is the mean value of the vector field at any space-time point (it should be kept in mind that p is invariant by 
time translation and so any mean value is stationary) is completely determined by the mean value at the origin of the 
coordinates, with the same density operator, modulo the replacement of thermodynamical parameters in (j28p . This 
particular value is strongly constrained by the symmetries of p. Let us identify the oj direction as that of the z (or 3) 
axis (see fig.H]) and consider the reflection with respect to z = plane and the rotation R3(7r) of an angle tt around 
the z axis; by repeating the same reasoning as for eq. (|T7l) for ^'^(0) we can easily conclude that the time component 
^"(0) is the only one having a non- vanishing mean value. Note, though, that the mean value on the right-hand side 
of (Uni) depends on the distance r from the axis because the density operator is modified by the replacement of the 
uniform temperature T with a radius-dependent Tq — 7T. Therefore, according to eq. (j29p and using (|25p. the mean 
value of the vector field can be written: 



V^{x) =ti[p:%{x)-] = ^(A-i)^otr 



-P°/To(r)+7(r)<.-J/To(r)+7(r)pQ/To(,.) . ^O(o) .1 = ^ y (3Q) 



i.e. it must be collinear with the four- velocity field u = (7, 7v) in eq. (^5)) and, therefore, its field lines are circles 
centered on the z axis and orthogonal to it. Similarly, we can obtain the general form of tensor fields of various rank 
and specific symmetry properties as a function of the basic four-velocity field. 

However, the previous derivation relies on the fact that the system is infinitely extended in space. Indeed, at a 
distance from the axis such that |o; x x| = 1 the velocity becomes equal to the speed of light and the system has 
a singularity. We cannot, therefore, take the strict thermodynamical limit V ^ 00 for a system with macroscopic 
angular momentum. Instead, we have to enforce a spatial cut-off at some distance and figure out how this refiects on 
the most general forms of vector and tensor fields. 

Enforcing a bounded region V for a thermodynamical system implies the replacement of all traces over the full set 
of states with a trace over a complete set of states \hv) of the fields for this region V, that we indicate with try: 

tr -J> try = '^{hv\ ■ ■ ■ \hv) 

hv 

The density operator is the same as in (jlSp with the partition function now obtained by tracing over the localized 
states. It may sometimes be convenient to introduce the projection operator: 



Py = J2\hv){hi 



which allows to maintain the trace over the full set of states of provided that we replace p with Pyp; for a generic 
operator A 

try [pi] = tr[Pypl] 
which amounts to state that the effective density operator is now pv'. 



Pv = ^Pv exp{-H/T + UJ-J/T + pQ/T) 

Zui 



(31) 



where: 



tr[Py eyiY>{-H/T + o; • J/T + pQ/T)] = try [exp(-i//T + • J/T + Q/T)] 



In order to maintain the same symmetry of the density operator in (jl8p . Py has to commute with J^, H, P^, the 
Lorentz boost along z and the reflection operator with respect to any plane parallel to z = 0, fl^ (see fig. [IJ. 
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These requirements are met if the region is a static longitudinally indefinite cylinder with finite radius R and axis 
uj, and we will henceforth take this assumption. 

There are two important consequences of having a finite radius R. Because of the presence of the projector Py , the 
previous derivation which led us to express the vector field according to the simple formula pop cannot be carried 
over to the case of finite (though macroscopic) radius. The reason is that Py does not commute with the Lorentz 
boost along v or, otherwise stated, a Lorentz boost along a direction other than z will not transform the set of states 
\hv) into themselves, as needed for completeness. So, one the one of the crucial steps in eq. (pS)) no longer holds and 
specifically: 



tr 



As a consequence, general vector and tensor fields will be more complicated than in the unphysical infinite radius case 
and get additional components. The most general expressions of mean value of fields in the cases of interest for the 
stress-energy and spin tensor will be systematically determined in the next section. The second consequence is that 
boundary conditions for the quantum fields must be specified at a finite radius value i?, but we will see that those 
conditions alone cannot ensure the validity of eqs. (ITT|) or (|12p . which are local conditions. 





n 



FIG. 1: Rotating cylinder with finite radius R at temperature T. 
the vectors of tetrad (|32p . 



Also shown the inertial frame axes and the spatial parts of 



V. TENSOR FIELDS IN AN AXISYMMETRIC SYSTEM 



In this section we will write down the most general forms of vector and tensor fields in an axisymmetric system, 
i.e. a system with the same symmetry features of the thermodynamical rotating system at equilibrium studied in the 
previous section. The goal of this section is to establish the conditions, if any, to be fulfilled by the superpotential to 
generate a good transformation of the stress-energy and spin tensors. 
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A. Vector field 

The decomposition of a vector field will serve as a paradigm for more complicated cases. The idea is to take a 
suitable tetrad of space-time dependent orthonormal four- vectors and decompose the vector field onto this basis. The 
tetrad we choose is dictated by the cylindrical symmetry: 

w=(7,7v) T—{"fv,^-v) n=(0,f) A: = (0,k) (32) 

where f is the radial versor in cylindrical coordinates, while k is the versor of the z axis, that is the axis of the cylinder 
(see fig. [I]). 

Due to symmetry for refiections with respect io z = const planes, the most general vector field V has vanishing 
component on k, and therefore: 

V = A{r)u + B{r)T + C{r)n (33) 

where A, B, C are scalar functions which can only depend on the radial coordinate r, owing to the cylindrical symmetry. 
Note the presence of two additional components with respect to the infinitely extended cylinder case in eq. pOI) . For 
symmetry reasons the only surviving component of the field at the axis is the time component, so i?(0) = C(0) =0. 
If the field is divergence-free, then C{r) = 0. 

B. Rank 2 antisymmetric tensor field 

Any antisymmetric tensor field of rank 2 can be decomposed first as: 

where: 

XP = - ie""^^ YP = AP'^u^ 

and, thus, X and Y are two space-like vector fields such that X ■ u — Y ■ u — Because of the reflection symmetry 
with respect to z = const planes, one has A^z = Ay^ = Aq^ = and this in turn entails that, being = 0, the only 
non-vanishing component of the pseudo-vector X is along k. Conversely, y is a polar vector and it has components 
along r and n which must vanish in r = 0. Altogether: 

AP'' = D{r)eP''P''kpU^ + E{r){rPu'' - t^uP) + F{r){nPu'' - n^uP) (34) 

with E{0) = F(0) = 0. Since: 

eP^P^kpU^ = nPT"" - n^rP 
(which can be easily checked), the expression ([M)) can be rewritten as: 

AP" = D{r){nPT'' ~ n^rP) + E{r){TPu'' - t^uP) + F{r){nPu'' - n^uP) (35) 

C. Rank 2 symmetric tensor field 

For the symmetric tensor SP'^ we will employ an iteration method in order to write down the most general decom- 
position. First, we project the tensor onto the u field: 

SP" = G{r)uPu'' + qPu" + q'^uP + QP" 

where q ■ u — and QP'^Ui, = 0. Then, we decompose the space-like polar vector field q according to ([55)1 

q = H(r)T + I(r)n 

with H{Q) ~ /(O) — 0, and we project the tensor Q in turn onto the vector field r: 

SP" = G{r)uPu'' + H{r){TPu'' + t'^uP) + I{r){nPu'' + n^uP) + J{r)TpT'' + hPr'' + /I'^r'' + S'^'^ 
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being h ■ u = h ■ t = (whence h = K{r)n with K{0) = 0) and ^^"^Ti, = 'B.^^Uy = 0. This procedure can be iterated 
projecting S onto n and the thus-obtained new symmetric tensor onto k. Thereby, we get: 

S^"" = G{r)u^u'' + H{r){T^'u'' + r^w'') + I{r){n^'u'' + n''u^') 

+J{r)T^T'' + K{r){n^'T'' + n'^r^) + L{r)n>'n'' + M{r)k^'k'' (36) 

However, since u'^u'^ — t^t'^ — n^n'^ — k^k" = g^"^ the last term can be replaced with a linear combination of all 
other diagonal terms plus a term in g^^'^ and the most general symmetric tensor can be rewritten, after a suitable 
redefinition of the scalar coefficients, as: 

S^"' = G{r)u''u'' + i/(r)(r^u'' + t^u") + /(r)(n^u'' + n'^u'') 
+J{r)T^'T'' + K{r){n''T'' + n^T^") + L{r)n''n'' - M{r)g'"' 

where H{0) = /(O) = K{0) = 0. 

D. Rank 3 spin-like tensor field 

The decomposition of a rank 3 tensor is carried out in an iterative way, similarly to what we have just done for the 
rank 2 symmetric tensor. First, we project the tensor onto the vector u and, taking the antisymmetry of fiiy indices 
into account, one obtains: 

^A,M^ = u^{f^'u'' - ru^") + u^V" + S^^u'' - E^'^u^ + T^'^"' (37) 

where all vector and tensor fields have vanishing contractions with u for any index. Particularly, using the general 
expressions (l33|) and ([35|) . the vector field / and the antisymmetric tensor F read: 

/ = £:(r)r^ + F(r)n^ ^ D{r){n^'T'' - n^T^") (38) 

with E{jS) = -F'(O) = 0. The tensor S can be decomposed as the sum of a symmetric and an antisymmetric part; 
having vanishing contractions with u, according to eqs. psp and (j36p . it can be written as: 

5]Am ^ N{r){n^T^' - n^r^) + P{r)T^T'' + Q(r)(n^T'" + n'^r^) + i?(r)n^n^ + S{r)k^k'' (39) 

wiht Q(0) = 0. The tensor T is projected in turn onto n and the above procedure is iterated. Then, similarly to 
eq. 1371: 

T^'^" = x'"'n^ + {h'^n'' - h''n^')n^ + Q^'^n" - G^'^n^ + A^'^"' (40) 

where all tensors have vanishing contractions with u and n. The antisymmetric tensor x must be orthogonal to u and 
n and, therefore, according to eq. psp . vanishes. On the other hand, the vector field h can only have non- vanishing 
component on t and so h = T[r)T with T(0) = 0. Finally, the tensor Q must be orthogonal to n, besides u, hence, 
using eqs. (|35p and dHS)) . can only be of the form: 

e^'^ = [/(r)r^r^ + V{r)k^k'' (41) 

Likewise, the tensor A can be decomposed onto r and, because of vanishing contractions with u and n, it can be 
written as: 

^ ^r{r)k^ (kf^r" - k^r^") (42) 

Putting together eqs. ((37)) . ((38)) . (p9)) . (|40)) . (|4T|) and p2|) . the general decomposition of a rank 3 tensor with 
antisymmetric iiv indices is obtained: 

$A,M^ ^ D{r){n''T'' - nV^)u^ + £;(r)(T^u'' - t^u^)u^ + F{r){n^'u'' - n''u^')u^ + iV(r)(n^T'' - 7iM)u^ 
-7V(r)(nV'' - n''T^)u'' + P{r)T^ {t^'u" - t'^u") + Q{r){n^T^' + n^^T^Yi" - Q{r){n^T'' + rfT^)u^' 
+R{r)n^{n''u" - n"u'') + S {r)k^ {k^" u" - k^u") + T(r)(T'^n'^ - T''n^')n^ 

+t/(r)r^(r^n'' - r'^n^') + V {r)k^ {k" n'' - k^n") + W{r)k^{k^'T'' - k^r") (43) 
with £^(0) = F(0) = Q(0) = T(0) = 0. 
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We are now in a position to find out the conditions to be fulfilled by the superpotential $ to be a good transformation 
of the stress-energy and spin tensors in a thermodynamically equilibrated system with angular momentum, as derived 
at the end of Sect. |IT1 

Let us start from eq. (jl2p . which is the most constraining. Since: 

we can write a rank 3 tensor ([43]) in the form of eq. ([T2|) as long as: 

V{r) = U{r) ^ F{r) 
P{r) = R{r) = S{r) 
E{r) = W{r) = -T{r) 

D{r) = N{r) = Q{r) = (44) 

which are definitely non-trivial conditions. If these are fulfilled, then the superpotential (|43p reduces to: 

$A,M^ ^ (i^(r)n'^ + E{r)T^' + P{r)u^')g^'' - {F{r)rf + E{r)T'' + P{r)u'')g^^' = K^'g^" - K^g^" 

Now, the field — F{r)n^ + E{r)T^ + P{r)u^ ought to be a constant one, according to eq. ([T3l) . Since its divergence 
vanishes, then F{r) ~ and, by using the definitions (|32|) . we readily obtain the conditions: 

F{r) = P(r)/7 = const E{r) = -P{r)ojr (45) 

In conclusion, only if a quantum superpotential is such that its mean value, calculated with the density operator 
([3T|l . fulfills conditions (|44|) and (f45|) . is the corresponding transformation ([3]) possible. Otherwise, the original and 
transformed stress-energy and spin tensors are inequivalent because they imply different values of mean energy, 
momentum or angular momentum densities. Since the most general form of the mean value of the superpotential, i.e. 
eq. (|43p is highly non-trivial, the inequivalence will occur far more often than equivalence. To demonstrate this, we 
will consider a specific instance involving the most familiar quantum field endowed with a spin tensor. 



VI. AN EXAMPLE: THE FREE DIRAC FIELD 



We now come to the possibly most significant result obtained in this work: the proof of a concrete instance of 
inequivalence, involving the simplest quantum field theory endowed with a spin tensor, namely the free Dirac field. 
It is well known that the from the lagrangian density: 

C = ^I'T^a^* - mf* (46) 

one obtains, by means of the Noether theorem, the canonical stress-energy and spin tensors p^ : 

^ i 

5A,M- = i¥{7\ 1:^"^ = ^^{7', [7^ 7l}* (47) 



where: 



crfc/2 

ak/2 



and (Tfe are Pauli matrices. The spin tensor obeys the equation 
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The couple of quantum tensors in (|47)) can be changed through the psuedo-gauge transformation in eq. ([3]). Ac- 
cordingly, if we take = 5, namely the superpotential as the original spin tensor itself, a symmetrized stress-energy 
tensor and a vanishing spin tensor are obtained: 



^ 7 

" 4 



- p_ - - ■• 



S'^^t"" = (48) 



This transformation is well known as Belinfante's symmetrization procedure. 

One may wonder whether these tensors, fulfilling continuity equations, still exist in a bounded region breaking the 
global translational and Lorentz symmetry, such as our cylinder with finite radius; or if, because of the boundary, they 
get additional terms with respect to the usual form. The problem of Dirac field with boundary has been tackled and 
solved by the authors of the MIT bag model [il| . First of all, it should be pointed out that the continuity equations 
^ certainly apply to tensors ([Tf]) and P5|) on-shell, i.e. for fields obeying the free Dirac equation within the cylinder. 
Furthermore, it is possible to find suitable boundary conditions, discussed in the next subsection, such that the fiuxes 
vanish, as needed, without introducing an ad-hoc discontinuity in the Dirac field. Thereby, the stress-energy and 
spin tensors retain the same form as in the usual no-boundary case and the integrals over the bounded region of the 
time components have the same physical meaning of conserved generators. It is also possible to derive field equations 
and canonical tensors (I47p from an action, what is shown in Appendix A. 

The spin tensor in eq. (|T7)) has a remarkable feature which makes it easier to check the equivalence of the two couples 
in eq. (|Tf|) and (|48p : because of the special properties of gamma matrices, the spin tensor is also antisymmetric in 
the first two indices: 

S^'t"" ^ (49) 

and thus the mean value of this tensor is greatly simplified. The antisymmetry in the indices (A, /i) dictates that all 
coefficients of symmetric A/i terms of the general form of this kind of tensor found in eq. (|43l) vanish: 

E{r) = F{r) = P{r) = Q{r) = R{r) = S{r) = T{r) = U{r) = V{r) = W{r) = 

and that D{r) — N{r), so that S is simply given by: 

S^^^"" = D{r)[{n^'T'' - n''Tf)u^ + [n^T^" - n'V^)!/,'^ - {n^r" - n''T^)u''] (50) 

and it is described by just one unknown radial function D{r). Therefore, according to the conditions (j44p . the 
Belinfante tensors are equivalent to the canonical ones (|T7)) only if D{r) — 0, i.e. only if the spin tensor has a 
vanishing mean value. 
For A = 0, eq. §Uii reads: 

and, because of the antisymmetry, the only non-vanishing components are those with both /i and v equal to 1,2,3, 
indices that we denote with We can then write, using ([5^ : 

= D{r)(v}V - v?v') = D{r)eijkk'' (51) 

Therefore, as expected, the time part of the spin tensor, contributing to the angular momentum density, is equivalent 
to a pseudo- vector field D(r) directed along z axis. 

According to eq. ([S]), the variation of energy- momentum density reads: 

^da ($"'°'' - - $'''"°) = ^da (5"^°'' - 5°^°'' - 5"'°'°) = -^daS^'""" (52) 

which implies at once that the energy density is unchanged because 5*^'"*^ = in view of (|49p . whereas the momentum 
density varies by a derivative. Using (j5ip and recalling the expression of curl in cylindrical coordinates: 

-'^ Belinfante ~ -"^canonical 2 "-'canonical ^ a ""^^ \ ' 

= ^c"anonical + 2^^'^^^)^ — ^c^anonical ~ 2 df^^^ ''^'^'^ 
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Note that this last equation imphes that the mean value of the canonical stress-energy tensor of the Dirac field has a 
non-trivial antisymmetric part if D'{r) ^ as, according to ([2]): 

Now we can write the angular momentum density variation by thermal-averaging eq. ^ with <f> = 5: 

•-^BclinfantG ~ "^canonical ^ 2 L ^ ^ ) X ^ )i 

The sum of all terms linear in the superpotential returns a — ^o./^" (gee eq. (jl2p ) while for the derivative terms we 
can use eq. ([5^ : 

= ^"a^nical " ^ ^5^5°'"'' - x'^d^S^'"^] S^'^'^ (54) 

Therefore, by plugging the expression of the mean value of the spin tensor in eq. (jSip . the angular momentum 
pseudo-vector corresponding to the angular momentum density in (|54p can be finally written: 

J'Bolinfanto = ^'canonical ~ ^ (^^ ^ ^^^dr^ ^ ^ ^ '^canonical - + -^(^)^ ^ (^^) 

In order for the canonical and Belinfante tensors to be equivalent, as has been mentioned and as it is apparent from 
eqs. ([55)) and ([55)) the function D{r) ought to vanish everywhere. If D'{r) ^ 0, the two stress-energy tensors give 
two different momentum densities and are thus inequivalent; if, on top of that, D'{r) ^ —2D{r)/r then the angular 
momentum densities are inequivalent as well. In the rest of this section we will prove that this is exactly the case, i.e. 
neither of these conditions is fulfilled. In order to show that this is not a problem arising from peculiar values of the 
field at the boundary, we will conservatively enforce boundary conditions such that the total energy, momentum and 
angular momentum operators obtained by integrating the fields within the cylinder are invariant by transformation 
([3]). Note that for this to be true, in the case under consideration, it is necessary that the function D{r) vanishes at 
the boundary, i.e. D{R) = 0, because the difference between total angular momenta is: 

f f f ^ dD{r) \ - f^°° d /r^ \ » 

J d^X (J'selinfante " ^'canonical) = - J dz d.^ dr r ( 2^~dr ^ ^^"^H ^ ^ J J '^^ 'dr \ ~2^^^'') ^ 

Thereby, we will demonstrate that, although the stress-energy and spin tensors in ()47p and ([^5)1 lead to the same 
quantum generators, their respective mean densities are inconsistent. The problem we are facing is then to solve the 
Dirac equation within a cylinder with finite radius and second-quantize the field. 



A. The Dirac field in a cylinder 



The problem of the Dirac field within a cylinder with finite radius has been tackled by several authors in the context 
of the MIT bag model [l2|. We first stress that, as we have done thus far, we take the viewpoint of an external inertial 
observer in Minkowski spacetime, seeing the spinning cylinder globally at rest. This observer can use either Cartesian 
coordinates or cylindrical coordinates to describe the system, the former being certainly more convenient to express 
tensor fields components while the latter are fit to solve the Dirac equation, as we will see. The most important issue 
in searching for a solution of this problem is the choice of appropriate boundary conditions, not an easy task because 
the Dirac equation is a first-order partial differential equation. The authors of the bag model [ll[ have shown that 
the following condition 

= m''7,,*(i?) = (56) 



^ Actually, in the paper [TH , the boundary condition chosen is = \I'(-R), but the change of sign is indeed immaterial. 
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ensures the vanishing of the fluxes in eq. ([6]) through the border and ahows non-trivial solutions of the Dirac equation 
within the cylinder (see fig.[T]) which, however, extend to the whole space without any discontinuity in the field. These 
boundary conditions, whence the whole solution, are not affected by the rotation of a possible material support defining 
the outer surface of the cylinder, as, for the inertial observer, the motion transforms the boundary into itself. The 
above equation, in the non-relativistic limit, entails the vanishing of the "large" components of the Dirac field at the 
boundary, that is one is left with the Schrodinger equation with Dirichlet boundary conditions. The eq. ([56]) implies 
that \E'(i?)^'(i?) = at the boundary [ll|, whence the vanishing of the outward current flux j^(_R)n^. Morevoer, since 

^'^'(i?) — 0, for any value of (p, z, t, the outer surface of the cylinder must be such that O^^'^'^Ir — S(i?) or: 

n^9^(¥^')(i?) = — = -S(i?) (57) 



•r 



r=R 



Thus, the flux of energy-momentum of the canonical tensor at the boundary is vanishing because, using eq. (|47p and 
eqs. ([56 1) .([ST ]) : 



JaV 2 Jgy 2 Jgy 2 J gy 2 J gy 

Likewise, for the orbital part of the angular momentum flux: 

dS x^'T^'^nx -{li^y) = ]: I dS x^(9'^(l^*) - {^l ^ u) = =^ f dS {x^'n" - x^n^') = 

dV 2 Jgy 2 Jgy 

where the last integral vanishes because of the geometrical symmetry z ^ —z. Finally, the flux of the spin tensor also 
vanishes at the boundary because, using (|T7| and ([5S)) : 

nxS^'^'^iR) = i (W^S^''* + *S^''^*) = (WS''^^' - = (58) 

Therefore, the eq. ([6]) applies and the integrals: 

= [ d^x f^" J^'' = [ d^x J°'^''' (59) 

Jv Jv 

are conserved. Since, we also have, from the Lagrangian, the usual anticommutation relations at equal times: 

{*,(t, X), ^l{t, X')} - SabS^i^ - x') {^a{t, X), *b(i, x')} = {^Ut, x), *^(t, x')} = 

it is easy to check that the conserved hamiltonian i/2 J d^x ^'^9*^ is indeed, as expected, the generator of time 
translations, i.e.: 

and, therefore, putting together the above equation with eqs. (|59)) and (|47)) we conclude that: 

[H, A] = (60) 

for the case under examination. 



The complete solution of the free Dirac equation for a massive particle in a longitudinally unlimited cylinder with 
finite transverse radius, with boundary conditions of the kind (|56p has been obtained by Bezerra de Mello et al in 
ref. and we summarize it here. In a longitudinally unlimited cylinder, but with finite transverse radius R, the field is 
expanded in terms of eigenfunctions of the longitudinal momentum, third component of angular momentum, transverse 
momentum and an additional "spin" quantum number p^ . The relevant quantum numbers n = (Pz, C(m,{,0' C) 
take on continuous (pz) and discrete values {M X{m i i) ■ The third component of the angular momentum M takes 
on all semi-integer values ±1/2, ±3/2, . . .; the "spin" quantum number ^ can be ±1 and the transverse momentum 
quantum number: 



C(M,i,i) = PTiR 



(61) 
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takes on discrete values which are zeroes, sorted in ascending order with the label I = 1,2,... and depending on M 
and of the equation: 



J\M-i\iPTR) + sgn(M)6W Ju,^i|(pTi?) = 



where J are Bessel functions and: 



m being the mass and: 



, (±) + S.'nT'T 
b\ = . 



(62) 
(63) 



^'(x) = ^C/„(a;)a„ + K.(x)&], 

n 

where On and 6n are destruction operators of quanta n while: 



the transverse mass we note in passing that ' = l/b[ . The Dirac field itself can be written as an expansion: 

(64) 



+ 00 



E-EE E / dp.^E EE/ 

n Af e=±lC(M.M) M ^=-1,11=1-^-^ 

The eigenspinors Un and Vn read, in the Dirac representation of the 7 matrices and in cylindrical coordinates {t, r, iy9, z): 

I •/|M_i|(PT/r) \ 



sgn(M)ACe6(+)j|^,+ i|(pTir)e'^^ 
HJ\M~\\^VTir) 

( J\M^^\^VTir) \ 

i sgn(M)K^&<-) J| 1 1 {pTir)e'^ 
KcJ|M+i|(PT;r) 
\ ~i sgn(M)6j-) J| 1 1 {pTir)e''^ J 



i[{M-l/2)^+p,z-et] 



2tt 



-i[(A/+l/2)(/)+p^z-et] 



2tt 



(65) 



with: 



and e = ^^p^ + + rn? being the energy. The eigenspinors (|65p are normalized so as to: 

/ d^x *t^ = Va],a„ + &„6t 
Jv 



that is with: 



d^x Ul{x)U^,{x) = / d^x V^{x)V^,{x) - 6„ 



d^x Ul{x)V^:{x) ^ 



(66) 



being (5nn' = Smm'^^^'^wS{pz — p'z) ^-nd the anticommutation relations of creation and destruction operators: 

{an, aj^/} = {6n, ^i/} = (^nn' {a„, 6„' } = {a]^, 6„'} ^ (67) 

The normalization coefhcient in (j65p obtained from the condition (j66p reads 13]: 



(C„) ' - 2TTR^J^^j_,^{pTiR)'-^^^{2R^mli + 2^MRmTi + mR) 



Pti 



(68) 



* In the rest of this section the symbol pTi stands for a discrete variable taking on (Af, ()-dependent values given by the eq. I|61|l or, 
later on, by eq. II79I I. 
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B. Proving the inequivalence 



For what we have seen so far, from a purely quantum field theoretical point of view, the Belinfante tensors (j48p 
for the Dirac field in the cylinder could be regarded as equivalent to the canonical ones in eq. (j47p because they give, 
once integrated, the same generators ([59l) . This happens because the condition ([7]) is met for $ = 5 (what follows 
from eq. ([55)) ) and this implies, taking eq. (|49l) into account, that all the integrands of ([7]) vanish at the boundary. 
Yet, these two set of tensors are thermodynamically inequivalent because, as it will be shown hereafter, it turns out 
that, using eq. (ISD) : 



S"''' = -try(p :vI/{70,S^^}*:) = D{r)e,,kk'' + ^ D(r) ^ 



(69) 



at some r ^ R (we have used the eq. ([ST]) ), with p written in eq. (IT51) . This will be enough to conclude that either 
the energy-momentum or the angular momentum densities or both have different values for different sets of quantum 
tensors, as previously discussed. Note that the boundary condition (j58p together with the general expression of the 
mean value of the spin tensor (j50p implies that D{R) — 0, i.e. its vanishing at the boundary. 
We can rewrite the inequality (1691) by taking advantage of the commutation relation: 



implying: 



S°^'^ = tTv[p ■■ :] - ig°hr[p : ^-7"* :] + ig"hT[p : ^'7'^* :] = try [p : ^^^'^^S :] ^ 



or, equivalently: 



D{r) = ^ezj35°''J' - ^tv[p -.^h.^sJ:'^-^ :] = tr[p :^ts3^:] ^ 



(70) 



where the indices i,j can only take on the value 1 or 2. In the above equation and henceforth, we can take the 
Heisenberg field operators at some fixed time t = because of the stationarity of density operator p. Hence we just 
need to show that: 



trv.[p:*^(0,x)I],^'(0,x):] ^0 



with 



/ 1 
0-100 
10 

Vo 0-1 



(71) 



(72) 



for some point x within the cylinder and our goal is achieved. 

To calculate the mean value of the spin density in eq. ([TT]) . we start by observing that (see Appendix B for the 
proof): 



tr[p 



n"n'J 



tv[pblbn'] = 



which allows us to work it out by plugging in there the field expansion ([64^ : 

1 



try[p:*t(0,x)E,*(0,x):] =^ 



Q{e-Mu+f,}/T ^ I 



[C/t(a;)S,[/„(x)] 



tiipa^bn'] = tr[p flnbn'] = (73) 
[V^{x)E,y^{x)]. (74) 



1 



where we have taken into account that the normal ordering of fermions is such that : b^bl^, := — fc|^/fon- By using 
eq. §B and ([72]) : 



An 
CI 
An 



i\{PTir) 



An « 
CI 
An 



M 

2 r 



i\{PTir) 



|M+i|(PTir)](l + 4) 



2 



^Ul{x)^,Un{x). 
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hence, by using eqs. (|68)) and (ffO)) . we can rewrite eq. ((74|) as: 
try[p:*t(o,x)S,^'(0,x):] =i?(r) 

1 1 



= EEE/ 

M {=±1 1=1 



Pti 



175) 



The mean value of the spin tensor is therefore given by the sum of a particle and an antiparticle term which are 
equal only for fi = 0. As expected, it vanishes for r = i? in view of the eq. (I62p . yet our goal is to show that it is 
non-vanishing at some point x not belonging to the boundary. It is worth pointing out that, if this is the case, the 
spin tensor has a macroscopic value because, as it is apparent from (j75p . it is proportional to the number density (in 
phase space) of quanta 1/ exp[(£ — Muj ± fi)/T + 1]. 

It is most convenient to consider a point belonging to the rotation axis, i.e. with radial coordinate r — because 
Bessel functions of all orders but vanish therein. By working out the eq. ((75|) . it can then be shown that D{0) ~ 
for uj = and that it is an increasing function oi uj/T thereafter (see Appendix C). It thence follows that 

^(O) ^ for uj/T > 

and therefore the inequality (|7ip must be true for small, yet finite, values oiui/T around the rotation axis. We note in 
passing that for w = the whole function D{r) must be vanishing because of symmetry reasons. In fact, if w = 0, the 
density operator pip enjoys an additional symmetry, that is the rotation of an angle tt around any axis orthogonal 
to the cylinder axis, say R2(7r). This transformation corresponds to flip over the cylinder, which leaves the system 
invariant provided that uj — 0, and has the consequence that any pseudo- vector field directed along the axis must 
vanish. 

Moreover, it is easy to show, again by using eq. (|75|) . that the derivative of the function D{r) vanishes in r = for 
it is proportional to terms, with > 0: 

2Jjv(0)j;r(0) = Jjv(O) (JAr-l(O) - JAr+l(0)) 

which all vanish because of the known properties of Bessel functions. Hence, the mean angular momentum density in 
r = differs between canonical and Belinfante tensors, i.e. rewriting the eq. (1551) for r ~ 0: 

J Bclinfantc(O) = canonical(O) — -D(0)k 

where £'(0) is finite for finite cj and positive. Thus, the Belinfante angular momentum density is lower than the 
canonical one by some finite and macroscopic amount. 

We point out that, had we used the definition of mean values ([1} without normal ordering, this conclusion would 
be unaffected. Indeed, the spin tensor (IT7|) is a bilinear in the fields and therefore the difference between the two 
definitions is a VEV of the spin tensor: 

tr(p :5:) = tr(p5) - (0|5|0) (76) 

Vectorial irreducible parts of the spin tensor, such as 5°'*^ , have a vanishing VEV if the vacuum is invariant under 
general rotations. The vacuum of the free Dirac field in the cylinder - defined by an|0) = - is indeed rotationally 
invariant. If degenerate vacua existed, the commutation of the hamiltonian with angular momentum operators that 
was shown before (see eq. (UHl)) would ensure that they belong to some irreducible representation of the SU(2) group. 
However, for the free Dirac field in the cylinder, the angular momentum operator along the z axis turns out to be: 

11 

SO that Jz|0) = on all possible degenerate vacua. This means that the only possible multiplet is one-dimensional 
and, thereby, the vacuum is non-degenerate and the second term in the eq. ((76|) vanishes. 



C. The non-relativistic limit 



It would be very interesting to calculate the function D(r) numerically to "see" the difference between the Belinfante 
and the canonical tensors and to make sure that this difference is not a rapidly oscillating function on a microscopic 
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scale, which would render the macroscopic observation of the difference impossible. This is, though, very hard in the 
fully relativistic case but relatively easy in the non-relativistic limit m/T 3> 1 because in this case the eq. yielding 
the quantized transverse momenta reduces to the vanishing of one single Bessel function. This happens because in 
the non-relativistic limit: 

,(+) imr+m \ IJlr+rn ^ 2_rn+^^ ^ 2rn ^ for ^ = 1 , , 

y^^> — ^ ^ — I PT PT PT ^ (77) 

« pj, I ni^^ ^ Pi/2m ^ ET^l for e = -1 



(78) 



so that the eq. ([62]) in fact reduces to: 



'j|,,^i|(pTi?) =sgn(M)^J|^_,|(pTi?) ^0 fore = 1 
J\ M- i I (ptR) = sgn(M) ^ J| , u{ptR)^Q for ^ = - 1 



Altogether, we can solve the equation Jl{ptR) = for all integers L — M + and take the quantized transverse 
momenta: 

PTi = %^ (79) 

where I — 1,2, . . . are now the familiar zeroes of the Bessel function of integer order L. 
We can now separate the particle and antiparticle terms in the eq. (j75p : 



Pti 



ii h J-^ e(^-*^"±'^)/^ + 1 8^2^J2 (pTii?)(2i?m|; + 2^MmTi + m) 



M C=±l 1=1 

with D{r) ~ D{r)^ + D{r)^ . In the non-relativistic limit one has: 

2Rm^ + 2^Mtot -f m ~ 2Rm^ + 2$Mto + m ~ 2Rm^ (80) 

where the last approximation is due to the obvious assumption Rm 3> 1 and that the term |^Mr7i| can be comparable 
to Rm? only if \M\ is very large. However, terms with large \M\ are either suppressed by the exponential exp[a;M/T] 
or by the Bessel functions, which effectively implements the semiclassical equality M « Rpx] since non-relativistically 
Rw? ^ RpTm Ri \M\m, the approximation (|80l) is justified. We then calculate the terms with = 1 ^ = — 1 
in the sum in eq. (|78|) separately. For ^ = 1 one sets Af + 1/2 = L and writes the integrand of eq. ([75]) . including 
approximation ()80|) and taking into account (|77p : 



1 p| JfL^MiPTir) - ^Jf,^{pTir) 1 1 J^ ip^.r) 



where pTi is a solution of the first equation in ((78)) . Similarly, for ^ = —1 one sets M — 1/2 = L and obtains, by using 
the second of the equations ([75]) : 

1 p't J\L\iPTir) - ||tj2^^^|(p^,r) 1 JfLiiPTir) 



g(e-L<^+c./2±M)/T + I IQn^R^m^ Jf.ApriR) Q{e^L.^-^/2±t.)/T _^ i IQ^2J12^2 2 ^ m 

(82) 



1 1 JfuiPTir 



Q{e-Lu,-u,/2±^,)/T ^_ 1 47,2^^2 j2^_^^| (p-^ji?) 

Now, by using approximations ([5n)) . (|8ip and we can write the non-relativistic limit of D{r)^ as: 

nirt = ^ f f r '^1^1^^"'^^ I ^1^1 (83) 

4^2^2 Z^y_ '^^-g(s-L^-^/2±M)/T^.l J2 (p^^^) g(.-L^+i./2±p)/T ^. 1 J2 (p^j^) ^^"^^ 

where the first term is to be associated to particles with spin projection -1-1/2 along the z axis and the second term 
to those with projection —1/2. Finally, the integral over pz in eq. ([55)) can be worked out by first introducing the 



19 



non-relativistic approximation e = m + p'^/2'm +p1/2m and then expanding the Fermi distribution. The final resuh 
is: 



L — — OO / — I Tl— 1 ^ 



J|i + l|(PTii?/?i) J|2-_^|(pT/i?/?i) j 

where we have purposely restored, for reasons to become clear shortly, the natural constants. 

It is very interesting to observe that the functions _D(r)^, hence D{r), are non- vanishing in the exact non-relativistic 
limit c — ?> CXI. Indeed, it can be seen from eq. (j84p that no factor ?i or c or powers thereof appear as proportionality 
constants in front of it, because the D{r)^ dimension is already - in natural units - that of an angular momentum; 
the only needed is in the exponent, which is compensated by a shift of the chemical potential, and the only fi's 
needed are those multiplying uj and in the argument of Bessel functions. Since h multiplies lu everywhere and D{r) 
vanishes for a; = 0, we also see that the difference between canonical and Belinfante densities is essentially a quantum 
effect, as it vanishes in the limit h 0; this is expected as the spin tensor exists only for quantum fields. 

For very small values of hu/KT these two functions are proportional to Tioj/KT itself since D{r)\^^Q — 0, as 
discussed in the previous subsection. Retaining only the n = 1 term of the series, corresponding to the Boltzmann limit 
of Fermi-Dirac statistics, and expanding the exponentials exp{±nhuj/2KT) at first order, one obtains the noteworthy 
equality: 

Dir)^ = ntripi-.^^i:^^:)^] ^ i|^to[p(:vl/tv,:)±] = n^-^ (84) 

where the superscript ± implies that one retains either the particle or the antiparticle term in the expansion of the 
free field and (dn/d^x)^ is, apparently, the particle or antiparticle density. The eq. (j84p can be shown by retracing 
all the steps of the calculations carried out for the spin tensor just replacing with the identity matrix. 

The function D(r) can be computed with available numerical routines finding a sufficient number of zeroes of Bessel 
functions, according to eq. (j78p . For the numerical computation to be accurate enough one has to make the series in 
L, I and n quickly convergent at any r. For the series in L, two requirements should be met: first (in natural units) 
u)/T <C 1 in order to keep the exponential exp[La;/r] relatively small and, secondly, the radius R should be such that 
Ry/mT is not too large; this condition stems from the fact that, as the Bessel functions effectively implement the 
semiclassical approximation \L\ ~ pxR and pr ~ V mT, the effective maximal value of L is of the order of RV mT. 
For the series in /, one has to set m/T ^ 1, so that large pr's are strongly suppressed; this is also the non-relativistic 
limit condition. For the series in n, one has to choose ^ so as to keep far from the degenerate Fermi gas case. The 
function D{r) as a function of r is shown in fig. [5] for ii = 0, R = 300, T = 0.01, to = 1 and two different values of w, 
10~^ and 2 • 10"''; the function {r/2)D'{r) — D{r), which is the difference between angular momentum densities for 
the canonical and Belinfante tensors, is shown in fig. |31 

The plots in figs. I2I3I show that the angular momentum density is larger in the canonical than in the Belinfante 
case almost everywhere, except for a narrow space near the boundary, whose thickness is plausibly determined by the 
microscopic scales of the problem (thermal wavelength or Compton wavelenght). Thereby, the observable macroscopic 
value of the differences between angular momentum densities, for a rotating system of free fermions, is the slowly 
varying positive one in the bulk. While the boundary conditions are needed to ensure the invariance of the total 
angular momentum, the rapid drop to zero within a microscopic distance from the cylinder surface tells us that 
the chosen boundary conditions at a mascroscopic scale of observation correspond to a discontintuity or a surface 
effect. Any macroscopic coaxial sub-cylinder of the full cylinder with a radius r < R will therefore have different 
total angular momenta whether one chooses the canonical or the Belinfante tensors in eqs. (l47l) and (l48t respectively. 
Such an ambiguity is physically unacceptable and can be solved only by admitting that these tensors are in fact 
inequivalent. 



VII. CONCLUSIONS AND OUTLOOK 



In conclusion, we have shown that, in general, couples of stress-energy and spin tensors related by a pseudo-gauge 
transformation ([3]) and allegedly equivalent in quantum field theory, are in fact thermodynamically inequivalent. The 
inequivalence shows up only for thermodynamical rotating systems and not for the systems - familiar in thermal 
field theory - locally at rest in an inertial frame. We have worked out exhaustively an instance of such inequivalence 
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FIG. 2: Function D{r), corresponding to the mean 
value of the canonical spin tensor for the free Dirac 
field, in a rotating cylinder at thermodynamical 
equilibrium as a function of radius r, in the non- 
relativistic limit. 



FIG. 3: Difference between the mean value of the 
canonical angular momentum density and the Be- 
linfante angular momentum density for the free 
Dirac field, in a rotating cylinder at thermodynam- 
ical equilibrium as a function of radius r, in the 
non-relativistic limit. 



involving the free Dirac field and shown that, surprisingly, the canonical and Belinfante tensors imply the same mean 
energy density but different mean densities of momentum and angular momentum. Particularly, the latter is almost 
everywhere larger in the canonical than in the Belinfante case for a small, yet macroscopic, amount. We would like 
to stress that this result does not to depend on an inappropriate treatment of the quantum field problem in a region 
with finite transverse size. First of all, field boundary conditions have been chosen so as to guarantee the invariance 
of the integrated quantities - i.e. the generators - and, secondly, the final spin tensor value ([M|) is proportional to 
particle density through fixed parameters uj and T. Hence, its non- vanishing does not apparently depend on spurious 
factors depending on the radius R which would probably be there if that was the result of an inaccurate calculation 
of finite size effects. 

What is the right couple of tensors? Needless to say, this is a very important issue; for instance, if it was found 
that the quantum spin tensor is not the trivial Belinfante one (i.e. vanishing) this would have major consequences in 
hydrodynamics and gravity, even more if its associated stress-energy tensor had a non-symmetric part, because this 
could imply a torsion of the spacetime (for a recent discussion see e.g. ref. [l3|)- Thus far in this work, no method has 
been discussed to answer the above question. Indeed, the problem of determining the right tensors can be approached 
both theoretically and experimentally. 

From an experimental viewpoint, in principle we could decide if a specific stress-energy or a spin tensor is wrong 
by measuring with sufficient accuracy the angular momentum density of a rotating system at full thermodynamical 
equilibrium kept at fixed temperature T and angular velocity uj. This measurement would, for instance, be able to 
reject the canonical or Belinfante tensor without even the need of resorting to relativistic systems as their difference 
has a non-vanishing non-relativistic limit, as has been discussed at the end of last section. In practice, at a glance, 
this measurement would not seem an easy one. According to eq. (|84p . in the non-relativistic limit the difference 
between these two tensors is of the order of hu/KT times h times the particle density, that is particles have a 
polarization of the order of huj/KT. This ratio is extremely small for ordinary macroscopic systems; assuming a 
large angular velocity uj, say 100 Hz, at room temperature T = 300 °K it turns out to be of the order of 10^^^. 
Notwithstanding, this is precisely the polarization responsible for the observed magneto-mechanical phenomena, the 
Barnett ^15] (magnetization induced by a rotation) and Einstein-De Haas (rotation induced by magnetization) effects. 
It is therefore possible that with some suitable experiment of this sort one can discriminate between spin tensors; this 
will be the subject of further investigation. We would like to point out that the effect could be enhanced lowering the 
temperature so much to increase the ratio hw/KT, e.g. with cold atom techniques. 

From a theoretical viewpoint, we cannot, for the present, determine a thermodynamically "best" couple of stress- 
energy and spin tensor. Yet, we can argue, on the basis of a thermodynamical argument - which can in principle be 
used to assess any other couple of tensors - that the canonical spin tensor is favoured over the Belinfante one; the 
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argument also elucidates why the source of the inequivalence is ultimately the second law of thermodynamics. Let us 
first write down the entropy of a system with cylindrical symmetry and large but finite radius R, using eq. (jlSp : 

S ^ -try[plogp] = logZ„ + i|l - |(Q> - ^ • (J) 

where ( ) stands for tr[p ]. All average quantities in the above equation, namely energy, charge and angular 
momentum, have a density, meaning that they are given by a volume integral of a supposedly objective function. Now 
let us assume that also entropy has a physically objective density, what is expected to happen for a macroscopic system; 
in this case, according to the above equation, the potential log Z^^ has an objective density as well. Furthermore, 
if the system is properly thermodynamical, the derivatives of this density dlogZ^/dF with respect to the intensive 
parameters 1/T,fi/T and uj/T should give the energy density, the charge density and the angular momentum density, 
similarly to what happens for the integrated quantities. 

Suppose that the system is initially non rotating, i.e. with u = 0. We have seen in Sect. IVII that in this case the 
function D{r) vanishes, hence there is no difference between the mean values of canonical and Belinfante tensors. Let 
us then focus on a coaxial sub-cylinder with radius r <^ R, yet large enough for it to be macroscopic. Keeping the 
temperature T and chemical potential ^ fixed, let us turn on very slowly a small angular velocity Aw. Thermody- 
namically, we can think of the sub-cylinder as "the system" and the rest as an angular momentum reservoir at full 
thermodynamical equilibrium with it as both have the same angular velocity. It is well known that a thermodynamical 
system at fixed T, V and /i maximizes the value of the thermodynamical potential log Z . Likewise, a thermodynamical 
system with fixed T,V,^ and ui maximizes the value of the potential \ogZ^. For the sub-cylinder with radius r, 
according to our previous thermodynamical assumptions, this potential is given by the integral over the region Vr of 
the density d\ogZ^/dV and its variation after switching on the rotational motion is: 

Vr,T,fl 

where {J)r is the angular momentum of the sub-cylinder. This variation should be maximal and since, for what we 
have seen in Sect. IVII in any properly macroscopic coaxial sub-cylinder with radius r < R, the canonical value of 
the angular momentum is larger than the Belinfante one, the final value of the thermodynamical potential of the 
sub-cylinder will be larger in the canonical than in the Belinfante case. Hence, if the system could choose between 
canonical and Belinfante tensors, the former would be certainly favoured. Of course, this method allows to discriminate 
between two couples of tensors but, for the present, does not permit to single out the "best" couple of tensors for a 
given quantum field theory. For instance, for the free Dirac field, it may happen that neither the canonical nor the 
Belinfante tensors maximize the thermodynamic potential of such a subsystem and another couple does. 



AlogZ, 



dlog Z^^r 

d{uj/T) 
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APPENDIX A - Action and canonical stress-energy tensor for the Dirac field in a bounded region 

Consider an action A of general fields V'" in a cylindrical region in fig. [1] and its variations: 



f f)r f)r \ j- f)r 



Particularly, if we take the Dirac action: 



A = I d^x ( ^*7''ap* - m^^j (86) 



and require it to be stationary with respect to variation of the fields with the boundary conditions (|56|) . we 

obtain the free Dirac equation. This can be shown by working out the derivatives in the boundary integral with C as 
in eq. (HS)) : 



n 



As the fields meet boundary condition (|56l) : 

as well as \I>^'(i?) = as a consequence of ([56|) . then 5'^'^ vanishes at the boundary and the second integral in eq. ((85|) 
vanishes. We are thus left with a bulk integral which has to vanish for general variations of the field, leading to free 
Dirac equations, as usual. 

For the calculation of canonical tensors from the action (|86|) . we follow ref. [16|] and use a space-dependent variation 
of the fields: 



(5*(x) = *(x + e{x)) ~ -^{x) ~ e^'{x)d^,^[x) 

with small e{x). This is a particular variation of the field which fulfills the condition (|56p if e{x) = at the boundary 
and this is what we set. If 4' is the solution of the equation of motion, then the variation of the action should vanish. 
After some easy calculations (see ref. [lij): 

The second term can be turned into a boundary integral which vanishes because e = there, as has been mentioned. 
The first term should also vanish and since e{x) is an arbitrary function, the divergence of the what can be easily 
recognized as the canonical stress-energy tensor must vanish. 

A similar reasoning leads to the conclusion that the angular momentum tensor is conserved. 

APPENDIX B - Calculation of the mean value of products of creation and destruction operators 

We follow the argument used in The aim is to calculate: 

try[pa]^a„/] 
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with p given by eq. (|18l) . For this purpose we define, with (3 = 1/T: 

aiil3) = e-P(^~^J-f^Q^aie^(^-^^-'^^^ (87) 
and similarly for an, bn and From the above equation it ensues: 

dam 



dp 



[al{^3),H-0JJ-^,Q] 



and, since: 



[H, al] = eal [J, a^] = Ma^ [Q, a^] = qa^ 
one readily obtains that eq. (|88p is equivalent to: 

which is solved by, being aJj(O) = aj^ : 

= aie-^'^'-^'^-'"'^ (89) 

We can now write: 

trvipalan'] = trv[pa],e^(^-"'^-'"5)e-^(^-"'^-^'3)a„,] = trv[e-^(^+-^'-^<5)a„, pale^'^^-^^-^Q^] 

where we have used the ciclicity of the trace, the definition of p in eq. and the eq ([57)) . It should be pointed out 
that the ciclicity of the trace can be used safely because a complete set of states for the cylinder with finite radius can 
be constructed with eigenvectors of the operators i?, Jz and Q; we could have also used the full trace and insert the 
operator Py discussed in Sect. HV] but this would have not changed the final result as this operator commutes with p. 
By using eq. (j89p and the anticommutation relation (j67p . the above equation can also be written as: 

try[pa],a„-] = try [pa„- a], (/?)] = try [pa„-a],]c-^(^-^^"-^«) = (-try [pa],a„-] + 5„„.) e-^(^-^^"-^«) 

whence: 

try[pa],a„,] = ^f^^.J^l^^^ ^ ^ 

The above method can be used for the calculation of other bilinear combinations of creation and destruction operators, 
leading to the equalities reported in eq. ([75]) . 



APPENDIX C - Calculation of the function D(r) on the rotation axis 



We calculate l?(0)|ij=o to show that it is vanishing as well as its derivative with respect to w/T to show that is 
strictly positive. The function D(r) in eq. (|75p is the sum of a particle D(r)^ and an antiparticle D(r)^ term: we 
focus on the particle term as the calculation for antiparticle is a trivial extension. Since Bessel functions of all orders 
but zero vanish (Jo(0) = 1), in r = in the sum of eq. ([75)) only terms with M — —1/2 and M — 1/2 survive: 



D(0) 



+ _ 



1 CXO „OQ 



^=±1 1=1 



j^g(e-i^+M)/T ^ Jo(p+,5i?^)2 [2r(pI ^ + m2) + ^^pl^^ + m^ + : 



+ m2 + m 



(90) 



24 



where we have defined p±.^ — ^'''^^''^ (see eq. (j6ip V We can rearrange the above sum by noting that the equation 
depending on indices (Af, ^) is the same for (— Af, — ^). In fact: 

J\-M-^\iC) + sgn(-Af)&W ^|-M+i|(C) - J\M+^\iC) - sgn(Af)&LV^|A/-i|(C), 

However, because of (|63p . —b'^:^^ = fo^ ■* = l/b''^\ and so multiplying the right hand side of above equation by 

sgn(M)&^^'' one gets the left hand side of eq. Hence, the zeroes of eq. (15^ and the one with "reflected" indices 

{—M, —^) are the same: 

C(-M-i,i) = C(M,i,i) (91) 

for any I — 1,2, . . .. Now we can redefine the indices in the second term of the sum in eq. (1901) by turning ^ into — ^, 
which changes nothing as ^ = — 1, +1 and write: 



D(0)+ 



1 pOO 



{=±1 i=i 



Jl(P-,-«i?)^ 



2i?( _^ + ) + £,\/p'- -f + ™^ + ^ 



We can replace P-,-^ with because of (|91l) and therefore: 
DiO)+ = 



8^ E E / dp J 

4=±i (=1 [ 



«2 



d2 



+ rn^ + TO 



(92) 



We are now going to prove that this latter expression is non-vanishing when w 7^ 0. First, we note that it does 
vanish when w = 0. In this case eq. (|92|) yields: 



D(0)+ 





Jl(p+^^i?)2-6^ 


V'jo(p+,4i?)2' 




[e(-+f)/T + 1] Ji(p+,^i?)2j„(p+^^i?)2 


2i? 


_t + to2 


+ i^p\,(^ +im? +171 



By using again (j9ip to replace with P-,-^ it is easy to show that the numerator of the integrand vanishes as: 



in view of the eq. (I62p . Therefore, the spin tensor density in r = vanishes for a non-rotating system, as expected. 
To show that it is no longer zero for a; 7^ we just need to show that the derivative with respect to uj /T in w = is 
not zero. One has: 



i+Y 



d 



d{u;/T) 



(e+p)/T 2 
6 P+^i 


JiiP+.^Rf- 








[e(-+f)/T + 1] Vi(p+,ei?)2Jo(p+,ei?)2 


2R 


P+A + 


+ W 


^ -I- to2 -|- To 



All terms are mainfestly positive except 
However, in this case: 



2R\ 



Pti 



£,^Jp^^j+rri? + rn in the denominator when ^ = — 1. 



2i?(p|; + TO,^) - + to2 + TO = ^p^i + to2 |^2i?yp|, + to2 - 1 j > yPri + C^Rm - 1) 

which is positive for a radius greater than half the Compton wavelength of the particle, that is positive for any 
actually macroscopic value of the radius R. The very same argument applies to the antiparticle term D(0)^ of the 
D{r) function in eq. (|75p with the immaterial replacement fi — >■ — /i, hence: 

^(0)U=o = --^DiO) >0 



dico/T) 



0=0 



